function [X,Y1,Y2,Y3]=demoData()
X0=1:40;
Y01=[0,1,3,2,4,5,7,6,21,15,...
    19,23,27,14,10,5,4,7,...
    5,8,3,9,11,22,29,31,...
    34,27,40,52,33,20,19,...
    16,14,60,55,54,66,70];
Y02=[35,34,38,32,26,17,25,8,15,2,...
    7,9,17,14,10,15,24,17,...
    25,28,33,49,31,42,29,31,...
    34,27,40,52,33,20,19,...
    26,24,10,15,14,6,0];

X=linspace(1,40,500);
Y1=interp1(X0,Y01,X,'cubic');
Y2=interp1(X0,Y02,X,'cubic');
Y3=X+3;
end